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Abstract 



We introduce a distinctive feature of spin-polarized transport, the Spin 
Coulomb Drag: there is an intrinsic source of friction for spin currents due 
to the Coulomb interaction between spin "up" and spin "down" electrons. 
We calculate the associated "spin transrestistivity" in a generalized random 
phase approximation and discuss its dependence on temperature, frequency, 
and electron density. We show that, in an appropriate range of parameters, 
such resistivity is measurable and propose an experiment to measure it. 
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Interest in spin-polarized transport has been growing dramatically in the last few years, 
spurred by the hope of realizing practical spin-electronic devices in a not too distant future 
In particular it has been shown that spin coherence can be maintained over large 
distances 8 S ~ 100/im and for long times T 2 ~ 10~ 9 — 10~ 8 s both in metals and in 
semiconductors ||. 

In this paper we introduce a distinctive feature of spin-polarized transport: in a conduc- 
tor, due to the Coulomb interaction, there is an intrinsic mechanism for friction between 
electrons of different spin, the "Spin Coulomb Drag" (SCD). For simplicity we shall re- 
strict our discussion to the case in which the spin state of each electron can be classified as 
"up" or "down" relative to the z axis. In the absence of impurities the total momentum 
P = Si Pi, where is the momentum of the i — th electron, is a conserved quantity. On the 
contrary, the "up" and the "down" components of the total momentum Pj = J2i PiT 1+ 2 Zi , 
and Pj_ = J2i Pij. 1 ~2 Zi i where a Z i is the the Pauli matrix for the z component of the i — th 
electron's spin, are not separately conserved even in the absence of impurities: Coulomb scat- 
tering can transfer momentum between spin up and spin down electrons thereby effectively 
introducing a "friction" for relative motion of the two spin components. If, for example, one 
of the two spin components is set into motion relative to the other, it will tend to drag the 
latter in the same direction. Or, if a finite spin current is set up through the application of 
an external field, then the Coulomb interaction will tend to equalize the net momenta of the 
two spin components, causing the difference (Pf) — (Pj) to decay to zero when the external 
field is turned off. 

The most dramatic manifestation of the SCD is the appearance of a finite trans-resistivity 
defined as the ratio of the gradient of the spin-down electro-chemical potential to the spin- 
up current density when the spin-down current is zero. This is completely analogous to 
the trans-resistivity measured in Coulomb drag experiments with electrons in two separate 
layers H-0], but in this case what makes the two electron populations distinguishable is 
not a physical separation but the different spin. In SCD the non conservation of the spin, 
caused mainly by the spin-orbit interaction, represents a "leakage" mechanism analogous to 
the interlayer tunneling in the usual Coulomb drag. 

First of all let us describe the SCD from a phenomenological point of view. Let Ej(t) 
and E|(t) be uniform effective electric fields that couple to spin j and spin { electrons 
respectively. These fields are sums of the ordinary electrostatic field plus the gradient of the 
local chemical potential, which can be spin- dependent. We restrict ourselves to the linear 
response regime, assume weak electron-electron and electron-impurity scattering, and ignore 
spin-flipping processes altogether. If v CT is the velocity of the center of mass of electrons of 
spin a, and N a the number of such electrons, the phenomenological equation of motion has 
the form 

mN a w a = -eN a E a + F ffS N a w a (1) 

where td is the Drude scattering time and F ct5 . is the Coulomb force exerted by spin <r(= — a) 
electrons on spin a electrons. By Newton's third law F ad . = — F Sa , and by Galilean invariance 
this force can only depend on the relative velocity of the two components. Hence, for weak 
Coulomb coupling we write 

n- 

F, Tff = -7miV ff -^(v ff -v ff ), (2) 
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where n a is the number density of electrons of spin a and n = n-f + nj_ is the total density. 
Eq. (HD defines the spin drag coefficient 7. Fourier transforming Eq. (HD with respect to 
time, and making use of the relationship j CT (c<j) = — en a v a (u) between current density and 
velocity, we obtain 

iujj a (uj) = + — 7 H ia(u) 

m \n t d J 

- -iyM- (3) 

n 

The resistivity matrix p a a i is defined as the coefficient of proportionality between the electric 
field and the current: = Pa,a'ja'- A quick comparison between this definition and 
Eq. @ shows that 7 is directly proportional to the spin trans-resistivity 

7 = -— Pn . (4) 

Let us now proceed to the microscopic calculation of the spin trans- resistivity. We start 
from the Kubo formula |7j for the uniform conductivity matrix 

a«A") = ---[nj.y + ((P<t;IV) H (5) 

iuj m m 

where ((A; B)) u represents, as usual f7], the retarded response function for the expectation 
value of A under the action of a field that couples linearly to B. The resistivity matrix is 
the inverse of the conductivity matrix. In the limit that both the electron-impurity and 
the electron-electron scattering are weak the P CT 's are almost constants of the motion and 
therefore ((P CT ; P<j'))w — ¥ 0. This means that the second term in the square bracket of 
Eq. can be treated as a small correction to the first |J. Inverting Eq. (|5|) to first order 
in ((P CT ; Pa'))uj and selecting the || matrix element we obtain 

= .% i^k. (6) 

e n ] n i 

It is convenient to recast this equation in a form that emphasizes the importance of the non 
conservation of Pj and Pj_. To this end we make use twice of the general equation of motion 

((A;B}}^ = -(([AB}) + z((A;B)U, (7) 

where A = i[A,H] is the time derivative of the operator A, and (..) denotes the thermal 
average. Thus, Eq. (|6|) can be rewritten as 

n (lA * ((P T ;P i )), + ^([P T ,P i ]) 

The commutator term controls the high frequency behavior of Pn(oo) and can be expressed 
in terms of ground-state properties 0. This term however gives a purely imaginary contri- 
bution to the trans-resistivity. Our present interest is in the real part of the trans-resistivity, 
which is controlled by the imaginary part of the force-force response function. 
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The force operator is given by 



i 2 

P - = ~y WqPmtP-q* VqP-q<r, ( 9 ) 

q q 

where v q = 47re 2 /g 2 is the Fourier transform of the Coulomb interaction, v e ~ % is the Fourier 
transform of the electron-impurity interaction, p qcr is the electronic spin density fluctuation 
operator, p q is the Fourier transform of the impurity density (a number), and V is the 
volume of the system. 

In calculating the force-force response function special attention must be paid to the 
contributions of the electron-impurity interaction. In the theory of the ordinary Coulomb 
drag JED] such contributions are zero on the average because the electrons in the two layers 
interact with two different sets of impurities, which are uncorrelated to each other. In the 
present case, however, electrons of opposite spin interact with the same set of impurities, so 
that electron-impurity terms, generated from the substitution of Eq. @ into Eq. (§), do not 
vanish upon disorder averaging. Happily, it turns out that these terms cancel out exactly at 
low frequency (u << Ep) and to leading order in the electron-electron and electron-impurity 
interactions [[II]]. Thus, the real part of the spin trans- resistivity takes the form 

Rep U (uj) = — — ^ — — VgVg, ■ 

n^n^uV 2 ^ 3 

J ™((P-qtPqi;Pq'TP-q'l))^ ( 10 ) 

where we have made use of the isotropy of the electron gas. 

We have calculated the four point response function 
X4 P (q, q',cij) = ((p-qtPqj.; Pq'tP-q'i))^ at finite temperature in a generalized Random Phase 
Approximation (RPA). The selected diagrams are shown in Fig. |l[ Because of its infinite 
range, the Coulomb interaction must be treated to infinite order, even when weak. The sum 
of the RPA diagrams has been evaluated by standard methods [[12] with the following result: 

1 n 2 0~< 8w - 1 1 

Re Pn (u, T) = L — £ q -v 2 ■ ^ L -l (11) 

n^n^Ve z q 3 y uj 

oo 7i \e~W -l)( e -/3(^') - 1) • 1 ] 

Here (3 = 1/ksT, with ks the Boltzmann constant, x"a'(lj u ) * s ^ e imaginary part of the 
RPA spin-resolved density-density response function, which is related to the noninteracting 
response function Xoo-(gS^) as follows 

bC l ((l^)} aa , = [X0cr] _1 (g,^)5crcr' ~ V q . (13) 

It is possible to show by simple but tedious algebraic calculations that this expression 
for the spin trans-resistivity pn(u,T) reduces, in the case of finite temperature and uj = 0, 
to the well know result of memory function and diagrammatic theories for the Coulomb 



drag [ rp | , |R|. Furthermore, for T = and tu 7^ 0, the RPA is equivalent to the decoupling 



approximation for the four-point response function used in |TJ| to calculate the dynamical 
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exchange-correlation kernel. Thus our calculation demonstrates that those two approxima- 
tions, quite different at a first sight, are simply RPAs performed in different limits. 

Let us focus on the low-temperature and low-frequency regime k B T « E F and 
uj « Ep, with Ep the Fermi energy. In this regime the imaginary part of the density- 
density response functions x'aa'iQi 00 ) * s a linear function of uj. In the limit of vanishing 
impurity concentration Xoa(q,w) is simply the Lindhard function, whose imaginary part, at 
low frequency, is given by Xoai^i u — ► 0) = — (m 2 /4ii)(uj/q) and whose real part can be 
approximated by its value at u — 0. Making use of this limiting form, the calculation of p^ 
can be carried in an essentially analytical fashion. The result is 



_ , T , ha 4ir 2 (k B T) 2 + uj 2 
Rep n {uj,T) - 



e 2 6(Ry) 2 
1 f°° dq 



247r 3 n i n T J-oo q 2 \e(q/a, 0)p 



(14) 



where a = h 2 /me 2 is the effective Bohr radius, Ry = e 2 /2a is the effective Rydberg, q = qa, 
n a = n a a 3 and e(q, uj) = 1 — v q Xo]{<l, uj) — VgXoii^ &) is the RPA dielectric function. Eq. ([14]) 
shows that, in the absence of impurities, p^(uj,T) is proportional to uj 2 for k B T << uj and 
to T 2 for uj « k B T. 

Modifications in the form of xoa(q,^) due to the presence of impurities can be taken 



into account through Mermin's approximation scheme |15[ . These modifications amount 
to replacing u/qvp by u/Dq 2 (D = v f t/3 being the diffusion constant) for uj < 1/r and 
q < 1/vpT where vp is the Fermi velocity and r is the electron-impurity mean scattering 
time. The uj and T dependencies of Eq. ( |i~4"D are not affected. 

Writing explicitly in Eq. (|TJ|) the dependence over r sa (where r sa = (47rn (T a 3 /3)~ 1 / 3 is 
the usual electron gas parameter) one can also see that p^{uj,T) ~ r st r s|5 so that p^ will 
strongly increase with decreasing electron density. In Fig. |2] we plot ph(uj = 0,T) for 
n-f = rii, at metallic densities (roughly 2 < r s < 6) in the temperature range 10/^ < T < 
70K. It can be seen that for temperatures of the order of 40 — 607^ (at which for example 
experiments on spin relaxation time using spin polarized currents have been performed 0), 
the spin trans- resistivity is appreciable (p-\i(oj — 0, T) ~ O.OlpVlcm). 

In the remaining part of this paper, we describe an experiment aimed at detecting the 
effect of the Spin Coulomb Drag and measuring the spin trans-resistivity. The setup is shown 
in Fig. 3: a paramagnetic metal film of thickness L is sandwiched between two ferromagnets 
polarized in the same direction. A battery is connected to the ferromagnets inducing a 
spin-polarized current from the first ferromagnet ("injector") through the paramagnet 
and toward the second ferromagnet ("receiver"). The injector and receiver are chosen to be 
semi-metals, i.e, they have only electron states of spin j at the Fermi level (see Fig.^[). It 
follows that the injected current j j is carried only by spin | electrons. If we choose L << 5 S , 
where 5 S the spin relaxation length, we can safely neglect spin-flipping processes and the 
polarized current entering the paramagnet will not relax before reaching the receiver. Spin 
relaxation lengths are relatively large in some materials (5 S ~ 100pm in Al 0]), so the 
condition L << 5 S is not particularly restrictive. Due to the SCD, the injected j| will drag 
spin | electrons toward the junction with the receiver. But, since there is no conduction 
band available in the receiver for spin j electrons the circuit will behave as an open circuit 
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for spin [ electrons, i.e., = 0. The vanishing of jj is an indication that the Coulomb drag 
force is exactly balanced by the gradient of the electro-chemical potential for spin down 

-eE i + m 7 ^ = 0. (15) 

where Ej = V/ij/e + E is the sum of the electrostatic field E and the gradient of the chemical 
potential What Eq. flTBD tells us is that due to the SCD there will be a measurable 
electro-chemical potential difference eE]l = erwyj^l/n^ for spin j electrons between two 
points within the metal separated by a distance / along the direction of the current. 

To measure this potential difference a second circuit including a voltmeter of very large 
resistance is connected to the regions of the paramagnet close to the junctions (See Fig. 3). 
Our purpose is to measure E^, so this second circuit must be driven by the spin j electro- 
chemical potential only. In order to accomplish this, we propose to use as contacts two 
semi-metallic ferromagnetic electrodes ("detectors"), similar to the injector and the receiver, 
but polarized in the opposite direction. In this way, for the same reasons explained before, 
the detection circuit will be "open" as far as spin f electrons are concerned, and the current 
flowing in the voltmeter will be exclusively driven by the electro-chemical potential difference 
of spin | electrons. The spin trans-resistivity will then be given by pfj_ = (AVd/Ii){A/1), 
where AVb is the voltage measured by the meter, A is the cross-section of the paramagnetic 
metal, / is the distance between the detectors, and if the current flowing between injector 
and receiver. As shown by our calculations, we expect a resistivity of the order of 10~ 2 fiflcm 
that is proportional to T 2 for k B T » uj. 

In summary we have pointed out a novel effect in spin-polarized transport, the Spin 
Coulomb Drag, and have proposed an experiment to observe it. We hope that this paper 
will stimulate experimental work aimed to the detection of this effect. 

This research was supported by NSF Grant No. DMR-9706788. We thank Shufeng 
Zhang for very valuable discussions. 
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FIGURES 



FIG. 1. The two series of "bubble" diagrams for the four-point response function \4p m the 
RPA. The vertices represent spin-density fluctuations p qa as labelled. 

FIG. 2. Density and temperature dependence of pn(0,T) in a paramagnetic metal. Top and 
bottom lines correspond to T = 70K and T = 10K respectively. Temperature is incremented in 
steps of 10K starting from the bottom. 

FIG. 3. (a) Experimental setup to detect the SCD effect: the voltage AV is applied between two 
parallel semi-metallic ferromagnets (injector (inj.) and receiver (rec.)) that sandwich a paramagnet 
(P). The voltage AVd is detected using two ferromagnetic electrodes (d) similar to the injector 
and the receiver, but polarized in the opposite direction, (b) Schematic bandstructure of injector, 
receiver, d and P. 
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Fig. 2 
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